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Abstract 

Lattice QCD simulations are now reaching a precision where isospin breaking effects 
become important. Previously, we have developed a program to systematically investigate 
the pattern of flavor symmetry beaking within QCD and successfully applied it to meson 
and baryon masses involving up, down and strange quarks. In this Letter we extend the 
calculations to QCD + QED and present our first results on isospin splittings in the pseu¬ 
doscalar meson and baryon octets. In particular, we obtain the nucleon mass difference 
of M„ - Mp = L35(18)(8) MeV and the electromagnetic contribution to the pion splitting 
- M„o = 4.60(20) MeV. Further we report first determination of the separation between 
strong and electromagnetic contributions in the MS scheme. 
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1 Introduction and general strategy 

Isospin breaking effects are crucial for the existence of our Universe. Our Universe would not 
exist in the present form if the n- p mass difference would only be slightly different. If it would 
be larger than the binding energy of the deuteron, no fusion would take place. If it would be a 
little smaller, all hydrogen would have been burned to helium. Isospin breaking in hadron masses 
has two sources, the mass difference of up and down quarks, and electromagnetic interactions. 
Both effects are of the same order of magnitude and cannot be separated unambiguously due 
to the nonperturbative nature of the strong interactions. This makes a direct calculation from 
QCD + QED necessary |[Ill2l[3]|. While substantial progress has been made, O is the only other 
published work to report simulations with fully dynamical QCD + QED. 

In dH [51 we have outlined a program to systematically investigate the pattern of flavor sym¬ 
metry breaking in three-flavor lattice QCD for Wilson-type fermions. Our strategy was to start 
from the SU(3) symmetric point with all three quark masses equal, m„ = = m^, and ex¬ 
trapolate towards the physical point keeping the average sea quark mass m = (m„ + + rus) /3 

constant. Eor this trajectory to reach the physical quark masses, m is tuned to the physical value 
of the average pseudoscalar meson mass = {m\o + j /3. We denote the 

distance from m by = niq - m {q = u, d, s) . This implies dm„ -l- dm^ -l- dm^ = 0 on our quark 
mass trajectory. To describe how physical quantities depend on the quark masses, we Taylor 
expand about the symmetric point [[51. This results in polynomials in m and dm^, which we clas¬ 
sify into representations of the SU(3) and S 3 flavor groups. As we keep m constant and change 
only the octet part of the mass matrix, to first order in dm^ flavor symmetry is broken by an SU(3) 
octet, leading to Gell-Mann-Okubo mass relations. We follow a similar approach here with QED 
added [| 2 l. 

The symmetry of the electromagnetic current is similar to the symmetry of the quark mass 
matrix. The simplifications that come from dmu + dm^+dnis = 0 in the mass case are analogous to 
the simplifications we get from the identity Cu + ed + es = 0. A difference between quark mass and 
electromagnetic expansions is that in the mass expansion we can have both odd and even powers 
of 6mq, whereas only even powers of the quark charges are allowed. We consider contributions 
of 0{eq) only. Hence, QED corrections can be simply read off from the mass expansion presented 
in dH, dropping the linear terms and changing masses to charges. 

Eor the masses of octet mesons with the flavor structure ah, and all annihilation diagrams 
turned off, we find to leading order in otem 

M^{ab) = Ml + a (dm^ -t 6mh) + (el + el +el) + jS™ (el + ej) + jd™ - e^f 
+ 7™ (el 5 mu elSnid + e^dmQ -t 7™ (e^dm^ -t e^dm/,) 

+ 7 ™ (ea - ebf (drua + dnib) + 7 ™ (el - ej) (dnia - drUb) 

+ 7™ ((^l + (^d + + dmb) -I- 7™ (ea + eb) (eudm,, + e^dtUd + e,dm,) 

up to corrections of 0(dml). Several of the coefficients in ([T]) can be matched up with different 
classes of Eeynman diagrams shown in Eig. [TJ The first diagram, with both ends of the photon 
attached to the same valence quark, contributes to (yS™ -l- >0™). The second diagram, with the 
photon crossing between the valence lines, contributes to yS™. The last diagram, with the photon 
being attached to the sea quarks, is an example of a diagram contributing to jS™. It would be 
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Figure 1: Examples of Feynman diagrams eontributing to the meson eleetromagnetie mass to 
Wavy lines are photons, eurly lines are gluons. 


missed out if the eleetromagnetie field was quenehed instead of dynamieal. Similar assignments 
hold for the mixed (eharge squared times mass) terms. For a single ehoiee of sea quark masses, 
the yS™ and 7 ™ terms ean be absorbed into the eonstant and the a term. However, for a 
eombined fit of both QCD and QCD + QED data we will need these eoeffieients. More details 
ean be found in Similarly, for oetet baryons with the flavor strueture aab we find to leading 
order in ctem 


M^iaab) = M, 

'EM ,2 , 2 , 2 ^ , oEM 


+ y 8 , 


0 + QTi ( 2 Sma + bnih) + 02 {Stria - Sttih) 

\el + el + e])+p\ 


{lei + el) + PT {ea - e,f + ySf - el) 


( 2 ) 


up to eorreetions of 0{5m^^. This exeludes the ease of baryons with three different quarks, as 
in the - A system 0. Again, the yS™ term ean be absorbed into the mass term M^. The 
eoeffieients /3™, yS™ and yS™ ean be matehed up with distinet elasses of Feynman diagrams 
similar to the ones in Fig. [B 

Our goal is to eompute the mass splittings of pseudosealar mesons and oetet baryons at the 
physieal point for QCD + QED. This amounts to determining the eoeffieients a, yS™ and 7 ™ 
in (HI) and dll)- It greatly helps to vary valenee and sea quark masses independently |151, whieh 
is referred to as partial quenehing (PQ). In this ease the sea quark masses remain eonstrained 
by m = eonstant, while the valenee quark masses yu^, yu^ and yu^ are uneonstrained. Defining 
duq = fig- in, the resulting modifieation of Eq. ([B to PQ oetet mesons is 

M^{ab) = Ml + a {6na + dyUy,) + yS™ {el + el + el) + /3™ {el + el) + /3™ (e„ - 
+ r™ {elSma + elSmd + + 7 ™ (e^dyu« + elS/Ub) 

+ r™ {ea - ebf {Sna + SiXb) + 7™ {el - el) (dyu^ - dyUy,) 

+ 7 ™ {el + el + el) (dyu„ + SjUb) + 7 ™ {ea + eb) + e,Sm,). 

For oetet baryons Eq. dB beeomes 


M^{aab) = Ml + ai {15na + d/iy,) + a 2 {5Ha - Sub) 


‘0 

oEM 

■Po 


/ 2 , 2 , 2\ , nEM /'-) 2 , 2\ , nEM / \2 , oEM / 2 2\ 

^s) + P\ + P2 - eb) + ^3 (e^ - e^) . 


(4) 


The eoeffieients a, P™ and 7 ™ in dl]) and dB are identieal to those in dB and d2l)- This is to say 
that hadron mass splittings are unaffeeted by PQ at this order, as PQ moves (e.g.) all oetet mesons 
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and baryons by the same amount. We would have to expand to eubie terms to see PQ errors in the 
splittings ill. Hence, PQ calculations offer a computationally cheaper way of obtaining them. 

In QCD + QED there is some ambiguity in the definition of the symmetric point. The defi¬ 
nition we have chosen is that the electrically neutral pseudoscalar mesons have the same masses, 
M^iuu) = M^(dd) = M^{ss) = M^{ds) = M^(sd) = M^(nn ), where n is a fictitious electrically 
neutral quark. As annihilation diagrams are neglected, different neutral mesons do not mix. We 
denote the Wilson hopping parameter k (introduced in (|7]) below) marking the symmetric point 
by Kg. We then have Strig = (trig - in) = 1/2/c®®^ - l/2k^ and S^g = (pg - in) = 1/2/c™' - \/2kg, 
setting the lattice spacing a = 1. It should be noted that even when all three quark masses are 
equal we do not have full SU(3) symmetry. Because of their different charges, the u quark is 
always distinguishable from the d and s quark. 

2 Lattice matters 

The action we are using is 

5 = 5g + 5a + 5“ +S^ + S^. (5) 

Here 5 g is the tree-level Symanzik improved SU(3) gauge action with gauge coupling jS = 6/g^, 
and Sa is the noncompact U(l) gauge action |[8ll3 of the photon, 

Sa = Yl Xi + H) - Ai,{x + v) - Ay{x)) . (6) 

We employ the nonperturbatively 0(a) improved SLiNC fermion action [fTOl for each quark 
flavor, 

= E E ~ l)e-‘’^‘'^>‘^^W,(x)q(x + fi) - q(x)(y, + l)e“^‘'^A^-f^Wl(x - fi)q(x - /))] 

+:^q(x)q(x) - ^Csw^q(x)o-i,yFgy(x)q(x)]. (7) 

^ fiv 

This action features single iterated mildly stout smeared QCD links with or = 0.1 |l5l and un¬ 
smeared QED links in the hopping terms, while the clover term contains unsmeared QCD links 
only. We keep the action deliberately local, as excessive smearing will lead to large autocorre¬ 
lation times. Stout smearing is analytic, so a derivative can be taken, which makes the HMC 
force well defined. The clover coefficient has been computed nonperturbatively in QCD IfTUII . 
We presently neglect electromagnetic modifications to the clover term. This will leave us with 
corrections of Oia^M a), which turn out to be no larger than the 0(a^) corrections from QCD 
in our simulations. We check this later by comparing neutral meson masses with different quark 
charges eg (Eig. [3]). Adding an electromagnetic clover term with c™ = 1 would leave us with 
corrections of DCctem g^a) (to this order in ctem), which is not a significant improvement, if at 
all. Simulations are performed using the HMC and RHMC IfTTlI algorithms. The gluon field and 
the EM field are updated sequentially. 

In this study, we limit our calculations to a single value of the strong coupling constant (lattice 
spacing) (5 = 5.50, where we have our largest sample of dynamical QCD configurations |[T2ll . 
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Furthermore, we restrict ourselves to simulations at the symmetric point, = 6md = 5ms = 0, 
which we define as X^JXl = 0.126, where = (mI + + M|_ + M|, + M|_ + M^) / 6 . We 

may use either X,^ or X^ to set the scale [|51. After several tuning runs carried out on 24^ x 48 
lattices we arrived at the k values Hu = 0.124362, Hd = Ks = 0.121713. At these k values, we 
study three different volumes, 24^x48, 32^x64 and 48^x96, with 0(2000) to 0(500) trajectories. 
We like to add that simulations at the symmetric point already catch the essential features of the 
physical QCD + QED vacuum, as flavor singlet quantities vary slowly along the in = constant 
trajectory flSl. 

On these ensembles we have computed PQ pseudoscalar meson and octet baryon masses for 
a variety of quark masses ranging from mps = 0.22 to 0.5, with Cq = -1/3,0 and +2/3. This 
leads to about 40 pseudoscalar masses and 70 baryon masses per ensemble. The baryons include 
several artificial states containing the fictitious n quark and charge 2 baryons with flavor structure 
uuu'. 

The action (|5]) is invariant under U(l) gauge transformations 

Af,(x) A^{x) + Af, a(x ), q{x) q{x). (8) 


However, this is not the case for propagators of charged particles, which demands fixing the 
gauge, as in perturbation theory. We choose the Landau gauge, which is defined by the condition 
AqA^{x) = 0, where A^ (A^) is the forward (backward) lattice derivative. The Landau gauge does 
not eliminate all gauge degrees of freedom, but allows for shifts A^a{x) of the photon field with 
A^a(x) = 0, where A^ = 0. To maintain (anti-)periodicity of the quark fields, a(x) must 

be periodic up to a transformation of the form 
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Ligure 2: Left panel: The background field on the 24^ x 48 lattice divided into three bins of 
approximately constant B^. Right panel: The bin averaged energy of the charged pion at rest 
against the bin averaged (B^). The black square (■) indicates the ensemble average (LnAv) of 
both aA/;r+ and (B^). The line is a one-parameter linear fit through the three mass points. 
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where is the extent of the lattice in /i direction. This gauge field redundancy can be eliminated 
by adding multiples of InjCqLf^ to A^{x), such that 
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<B^< 


7T 

\eq\Lf, ’ 


B,= 


1 

V 




(9) 


Taking = -1/3 in ([91) serves both charges. The advantage of this procedure is that it leaves 
the fermion determinant and Polyakov loops for all quark flavors unchanged. In other popular 
gauges |l3l this is not the case, but results in a permanent Polyakov loop 77^^ L™(Xy) = 1, 
which we would not know how to correct for in a simple manner. 

The constant background field can be factored out from the link matrices, configuration by 
configuration, and absorbed into the quark momenta by straightforward algebraQ This leaves us 
with photon propagators that are devoid of zero modes. In the presence of a constant background 
field 5^ the correlator of a single hadron H thus becomes ||9ll 




( 10 ) 


where M^, p and eu are mass, three-momentum and electric charge of the hadron, respectively. 
This amounts to a shift of the rest energy of the charged hadrons, Mh ^ Mh + 

e\ B^/2 Mh. We determine the ensemble average of directly on each of our three volumes. The 
result is 0.024, 0.0079 and 0.000095 for the smallest to largest volumes, respectively. To extract 
masses, we remove the influence of the background field effect by subtracting the associated 
kinetic energy from the ensemble averaged lattice energy. To demonstrate the validity of this 
procedure, we have divided a subset of our 24^ x 48 ensemble into three bins of approximately 
constant background field in Fig. [2] and plot the corresponding lattice energies for each of these 
bins against the corresponding B^. It shows that both the energies of the individual bins as well as 
the ensemble averaged energy fall on a single straight line, in line with our subtraction method. 
On the 48^ x 96 lattice the effect of the background field is comparable to our statistical precision. 
With the zero modes removed, we then can employ established methods, such as [fTSll . to correct 
for the remaining electromagnetic finite size effects associated with the long-range tail of the 
photon field. Any residual effect of the background field will only act to modify the recoil energy 
of any charged hadron propagator within loops. 

Our strategy is to simulate at an artificial coupling = 1.25, and then interpolate between this 
point and pure QCD to the physical fine structure constant ctem = 1/137. This value is chosen so 
that electromagnetic effects can be easily seen, but is still small enough that they scale linearly in 

and we do not need to consider higher order terms. Most importantly, Z 3 = 0.94(3), obtained 
from the vacuum polarization. Furthermore, in Figs. 4 and 5 of |l 6 l we have plotted l//c^ , l/Kq 
and the bare quark mass at the symmetric point, 1 jlRq-l jlK^, against for Cq = -1/3,0 and 2/3 
and found that all three quantities lie on a straight line. In addition, we find that the coefficients 
a and a\, a 2 m ^ and dH) agree to a good precision with the corresponding numbers in pure 
QCD [l5]|. This rules out significant higher order corrections in e^. 

'Note that the transformation q{x) —> q{x), q(x) —> q(x)e~'^^ amounts to a shift of quark momenta p —» 

p +A. 
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Figure 3: The neutral pseudoscalar meson masses Mpg(aa) on the 48^ x 96 lattice as a function 
of d/if for quark charges = -1/3 (blue), 0 (green) and 2/3 (red). Solid (open) symbols refer 
to unitary (PQ) masses. The horizontal lines display the physical n^,K^ and tJs meson masses, 
the vertical lines indicate the physical u and s quark masses. 

3 Results 

After the initial small volume tuning runs, it turns out that the chosen k values do not quite 
satisfy our constraint of equal neutral pseudoscalar meson masses. A more accurate estimate can 
be determined from a fit to the pseudoscalar meson masses. On the 48^ x 96 lattice we obtain 

k„ = 0.124382, = = 0.121703, k„ = 0.120814, (11) 

which is only a small displacement from the underlying simulation kappas. We shall expand 
about these k values in our subsequent fits. 

In contrast to QCD, equal meson masses at the symmetric point no longer mean equal bare 
quark masses. We renormalize the quark masses to remove this defect. We do so by absorbing the 
QED terms of the neutral pseudoscalar mesons into the quark self-energies. On our symmetric 
background, = dm^ = dm^ = 0, this is achieved by replacing d/Uq by the ‘Dashen’ scheme 
mass [|6l 

d//f = [l+(r™/a)e;]d//,. (12) 

Substituting (fT^ into ([3]), and absorbing into and y™ into a, we obtain in the ‘Dashen’ 
scheme 

M\ab) = + a (d/if + d/if) + (e, - e,f 

+ r™ (e. - e,f (6iu^ + + r™ (e^ - el) (d/if - d/if). 

Note that since we choose the neutral pseudoscalar mesons to have the same mass, = 0 by 
definition. We define the critical point, for each flavor, to be the point where the masses of 
the neutral pseudoscalar mesons vanish. It is then easily seen that the ‘Dashen’ scheme quark 
masses are all equal at the symmetric point, /if = q = u,d,s and n, see O for further 

details. It follows that the total electromagnetic contributions to the neutral pseudoscalar meson 
masses, M„o and M/^o, are zero. In Fig. [3] we show the neutral meson masses M^{ad) against 
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Figure 4: The eharge (p and baryon masses Mjj(aab) on the 48^ x 96 lattiee as a 

funetion of 26p^ + 5p^. Solid (open) symbols refer to unitary (PQ) masses. 

d/if. It is striking to see that the data fall perfectly on a straight line, which strongly supports 
our group-theoretical classification of SU(3) flavor symmetry breaking as well as octet (Gell- 
Mann-Okubo) type mass splitting. 6 >(q;em a) corrections would result in deviations from the 
straight line proportional to ef for which we see no evidence. To be consistent, we also expand 
the baryon masses in terms of the ‘Dashen’ masses, 

M\aab) = Ml + (2d/if + d/if) + ^2 (d/if - dyuf) 

+ /S™ (2i.2 + el) + {ea -e,f+ {el -el). 

In Fig. m we show the charge e^ = +\ baryon masses M^{aab) against 2d/if -l- d/if. Again, the 
data fall perfectly on a straight line, in accord with our flavor expansions. 

For the total contribution of QCD -l- QED it does not matter which scheme we use to define 
the quark masses, but for the individual contributions of QCD and QED it will make a difference. 
The fits of (fTSl) and (fT4l) to the lattice data are quite robust, giving;^^/dof = 0.7 - 1.2. To obtain 
physical numbers, we extrapolate the coefficients yS™ and y™ to ctem = 1/137 by scaling them 
with a factor ~ 10/137. In our extrapolation to the physical point we keep the sum of the quark 
masses constant. We choose and -l- to determine the physical k values. 

In Fig. [5] we show the result of the fit to the meson and baryon masses on the 48^ x 96 lattice. 
We obtain = 0.128(3), which is to be compared with the physical value, 0.126. This tells 

us that we have hit the symmetric point with remarkable precision. Using X„ to set the scale, the 
lattice spacing turns out to be a = 0.068(2) fm. The figure also indicates that the baryon masses 
extrapolate nicely to their experimental values, leaving little room for quadratic terms. Similarly 
good results are found on the 32^ x 64 lattice. Having found the k values of the physical point 
and the point where the ‘Dashen’ scheme masses vanish (the critical point), we can determine the 
quark masses. For the quark mass ratios we find on the 48^ x 96 lattice in the ‘Dashen’ scheme 

— = 0.52(2), — = 19.7(9). (15) 

rud nid 

In we have shown how to switch between the ‘Dashen’ and MS schemes. Applying this, we 
find the ratio in the MS scheme at //^ = 4 GeV^ decreases by less than a percent, whereas 
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Figure 5: Fan plots of pseudoscalar meson (top) and baryon masses (bottom) on the 48^ x 96 
lattice as a function of -l- d/ij, with -l- d/ij + S/Us = 0. The baryon masses are the averages 
of the isospin doublets. 


nijm-d remains a renormalization group invariant, even in the presence of QED. Hence Eq. (fTSl) 
represents our results in the MS scheme at f/ = 4 GeV^. 

In this Letter we are primarily interested in the isospin splittings of pseudoscalar meson and 
octet baryon masses. To get to our final numbers, we need to correct for finite size effects first. 
From QED we expect power-law corrections, due to the photon being massless, in addition to 
exponential corrections from QCD. We correct for QCD finite size effects by using the results 
of lfT4l[T5]l . adapted to three flavors of PQ quarks. In case of the nucleon the corrections amount 
to approximately 1% on the 48^ x96 lattice and to 5% on the 32^ x64 lattice. Having successfully 
removed the zero modes, we can correct for the remaining QED effects by employing the mass 
shift formulae of effective field theory (EET) lfT3ll . We test this in Eig. 0 where we compare 
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Eigure 6: The QED contribution to the p - n mass splitting on the 32^ x 64 and 48^ x 96 lattices 
compared with the prediction of lfT3ll . 
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Figure 7: Mass splittings of pseudoscalar meson (top) and baryon masses (bottom) as a function 
of the spatial size of the lattice. The numbers on the largest volume have been extrapolated to 
infinite volume using lfT3]l . 


the QED contribution to the p - n mass splitting, {Mp - M„)qed, with the prediction of [[T3ll on 
our two largest volumes. We find good agreement between the data and the analytic expression, 
indicating that QED finite size effects are well accounted for by EFT. In Fig. |7] we present our 
QCD + QED results for the isospin splittings of mesons and baryons as a function of lattice size. 
The curves represent the predictions of lfT3l . They have been drawn through the points on the 
48^ X 96 lattice. We find good agreement between the curves and our points on the two largest 
lattices, while the data on the 24^ x 48 lattice (with L « 1.6 fm) appear to lie outside the range of 
validity of the expansion. We consider the extrapolation of the 48^ x 96 lattice points to a/L = 0 
by lfT3]l our best estimate of the infinite volume result. We compare this result with a fit to the 
points on the two largest lattices. The differences are taken as an estimate of systematic error. 
In Table [U we list our final results for the mass splittings in the infinite volume, for the total and 
the QED contribution separately. Following |I3, we find the QED contributions in the ‘Dashen’ 
scheme and the MS scheme at = 4 GeV^ to differ by less than a percent. As a result, the 
QED contributions in Table [U also represent our results in the MS scheme at = 4GeV^. The 
traditional way of expressing the electromagnetic contributions is through and 

the e parameter, 

(M|, - m2„)qed - Ml. + = eA, . (16) 

On the 48^ x 96 lattice we find 6 = 0.49(5), which translated to MS gives ll6ll 

6 = 0.50(6). (17) 

This result is well within the range quoted by FEAG ifTTll . albeit with significantly reduced 
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AM 

QCD + QED 

QED 

QCD [m 

Experiment 

M„+ — Mj^u 


4.60(20) 


4.59 


4.09(10) 

-1.66(6) 


3.93 

Mn-Mp 

1.35(18)(8) 

-2.20(28)(10) 

3.51(31) 

1.30 

M^- — M^+ 

7.60(73)(8) 

-0.63(8)(6) 

9.07(47) 

8.08 

Ms - Mso 

6.10(55)(45) 

1.26(16)(13) 

5.58(31) 

6.85 


Table 1: Mass splittings in the infinite volume, in units of MeV. The first error is the statistical 
error from the extrapolation of the points on the 48^ x 96 lattice. The second error (if any) 
is a systematic error estimated from the fit to both the 48^ x 96 and 32^ x 64 volumes. The 
QCD + QED and QED results are compared with previous results from pure QCD [fT^ and the 
experimental numbers. 


uncertainty. We now can compare the baryon mass splittings of this calculation with our recent 
results from pure QCD [fT^ . The QCD numbers are quoted in the fourth column of Table [B 
They have been brought in line with our new value of e (fT71) . Both sets of results are found 
to be largely consistent. It is worth emphasizing that the QED and pure QCD contributions 
to the nucleon mass splitting sum up nicely to the total QCD + QED contribution, which is 
encouraging. Einally, in the last column of Table [B we quote the experimental mass splittings. 
We observe good agreement for both octet pseudoscalar mesons and octet baryons. Since we 
have not yet computed the QCD contribution to the mass from n^-ri mixing, arising from 
quark-line disconnected diagrams, we only quote the QED contribution to the M„+ - mass 
difference. It is worth noting that phenomenological estimates for the disconnected contribution 
are of the order of 0.1 MeV lUSll , which is within the precision of our present calculation. EigureH] 
summarizes our results. 

Both the total QCD -l- QED mass splittings as well as the QED contributions satisfy the 
Coleman-Glashow relation [[T9ll by construction. So do the experimental values, which once 
again supports our group-theoretical approach and truncation (fT4l) . The QED contribution to 



Eigure 8: Mass splittings AM of octet pseudoscalar meson and baryon masses compared to 
experiment. 
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Figure 9: The allowed ratio of quark masses niulnid for a range of against qtem- The solid eirele is 
our result in the MS seheme. The region of no fusion is to the left, the region where all hydrogen 
is converted to helium stars is to the right. 

the n - p mass splitting in the ‘Dashen’ and MS schemes turns out to be somewhat larger (in 
absolute terms) than the numbers derived from the Cottingham formula [|20l . It should be noted 
though that the individual estimates GOl cover a wide range of values. To accommodate the 
lower numbers from the Cottingham formula, the result of pure QCD ifT^ (fourth column of 
Table [B would have to be smaller by a factor up to two as well. Our QED result is also larger 
than the recently reported lattice number in |I21 . In our approach the QED and QCD separation 
is defined within the meson sector. In contrast, [[3]| chose the QED part of the - Z~ mass 
difference to be zero, for which we identify a clear nonzero signal. This would be the case 
if (2/3)jS™ + yS™ + (l/3)/3™ = 0 in our mass expansion (fT4l) . A fit to our data with this 
constraint gives (M„ - Mp)qed = -1.71(28)(10) MeV in the ‘Dashen’ scheme. While this result 
is largely compatible with the analysis of Walker-Eoud, Carlson and Miller If20l . (M„ -Mp)qed = 
-1.30(50) MeV, it illustrates quite clearly that the QED part of the n-p mass difference depends 
sensitively on how electromagnetic and strong contributions are separated. While our results do 
not support higher order terms in the quark mass expansion, it may be possible that one source 
of the discrepancy could be related to nonlinearities in the chiral behavior of the electromagnetic 
self energy ll2TI that are not being captured by the Taylor expansion. 

As discussed in the introduction, the existence of the Universe as we know it is highly sensi¬ 
tive to the magnitude of the n-p mass difference. Having an analytic expression for the mass 
of neutron and proton, Eq. (fT4l) . we can express the allowed region in terms of the fundamental 
parameters m„, and aEM, as shown in Eig. [9l Not shown are the bounds on ctem from the 
stability of atoms If22]| . It turns out that both ctem and the ratio of light quark masses m„/m^ 
are finely tuned. At the physical fine structure constant the ratio is restricted to a narrow region 
around m^/m^ = 0.5. 
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4 Conclusion and outlook 


We have outlined a program to systematieally investigate the flavor strueture of hadrons in a full 
QCD + QED lattiee simulation. By treating the valenee quark masses differently to those in 
the sea allows for a range of valenee quark masses and eharges to be explored and signifieantly 
enhanees our ability to aeeurately eonstrain the fit parameters in our flavor-breaking expansions. 
As a result, we have sueeessfully eomputed the isospin splittings of pseudosealar meson and oetet 
baryon masses. By using our reeently introdueed ‘Dashen’ seheme as an intermediate step |j 6 l, 
we are able to quote the first lattiee results for the QED eontribution to the n - p mass splitting 
in the MS seheme. 

The ealeulations have been done at lattiee spaeing a = 0.068 fm. At this lattiee spaeing 
diseretization errors are expeeted to be less than 2 % IfT^ . whieh are well below our present 
statistieal and systematie errors. To reduee the errors and gain full eontrol over the infinite 
volume extrapolation, simulations on 64^ x 128 lattiees and larger will have to be done. To 
further eonstrain our fits, and test for potential effeets, we have started dynamieal 1 -l- 1 -l- 1 
flavor simulations along the m = eonst line, with dm„ dm^ 0 and sea quark masses 

approaehing the physieal point. Einally, future simulations will also naturally be required on 
lattiees with different lattiee spaeings to allow for a eontinuum extrapolation. 
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